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Abstract: The coordinate-space Hartree-Fock-Bogoliubov (HFB) approach with quasiparticle blocking has been 
applied to study the odd-A weakly bound nuclei ^^Mg, in which halo structures have been reported in 

experiments. The Skyrme nuclear forces SLy4 and UNEDFl have been adopted in our calculations. The results with 
and without blocking have been compared to demonstrate the emergence of deformed halo structures due to blocking 
effects. In our calculations, ^®B and ®^Mg have remarkable features of deformed halos. 
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1 Introduction 

Nuclei with extremely unbalanced neutron to proton 
ratios, close to drip lines, can have a lot of exotic phe¬ 
nomena such as halo states and pygmy resonances[T]. 

From the theoretical point of view, drip-line nuclei are 
weakly bound superfluid quantum many-body systems 
and can offer a unique environment to test threshold ef¬ 
fects as well as to investigate effective interactions. In 
the experimental aspect, the new-generation radioactive- 
isotope beam facilities such as ERIE [5] and HIAF [3] are 
being built to explore new physics near drip lines and to 
study nuclear processes in astrophysical environments. 

One of the most important reasons to study ex¬ 
tremely neutron-rich nuclei is their exotic surface prop¬ 
erties and dynamics, due to the intertwining of the halo 
decoupling effect, surface deformations, pairing correla¬ 
tions and continuum couplings [IHS]. In particular, the 
deformed halo structure can have decoupled core and 
halo deformations and various situations can happen [7] . 

In a series of previous studies based on the coordinate- 
space Hartree-Fock-Bogoliubov approach, we have stud¬ 
ied exotic deformed halos and collective modes in even- 
even weakly bound nuclei [HE]. Actually deformed halos 
in experiments up to now are all odd-mass nuclei such as 
^^Be |8], ®^Ne [9] and ®’^Mg [10]. Theoretical studies of 
the odd-A nuclei is not so popular because of additional 
computational complexity when dealing with pairing cor¬ 
relations. Indeed, the quasiparticle blocking effect has to 
be taken into account for spin-polarized systems such as 
nuclei and cold atomic gases m- 

In the past few years, several odd-N halo nuclei with 
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deformations have been reported experimentally, such 
as ®^Ne [9] and ^^Mg [TO], implying more opportunities 
than expected. Deformations and weak-binding effects 
together can change the standard Nilsson shell evolu¬ 
tion and drive s or p orbits that are not even around 
Fermi surfaces into halo states. Very recently, the de¬ 
formed halo ®^Mg was observed experimentally and cur¬ 
rently ^^Mg is the heaviest halo nucleus [10]. A large 
p-wave component is required to reproduce the experi¬ 
ment. Theoretical studies are few since the experiment 
result of ^’’Mg has come out only very recently. In addi¬ 
tion to the odd-N halo nuclei, we are also interested in 
the odd-Z halo nuclei such as ^^B and «B[I2]. In fact 
the neutron halo could also be influenced by the polariza¬ 
tion effect due to the blocking of proton orbits. Several 
experimental and theoretical studies suggest ^’^B has a 
halo structure |13H18j . For ^®B, the experimental stud¬ 
ies suggest it may have a 4n-halo structure [12]. Recent 
experiments found the valence neutron of ^®B has large 
d-wave component which hinders the formation of the 
halo structure [13]. Therefore the purpose of this work is 
to employ the Hartree-Fock-Bogoliubov(HFB) approach 
with quasiparticle blocking in deformed coordinate space 
to study these odd-A halo nuclei. 

HFB theory can properly treat the coupling be¬ 
tween the weak-bound orbits and continuum spectrum, 
which is essential for weakly bound nuclei. The con¬ 
tinuum states affect not only the properties of ground 
states but also excited states [5] [20]. Solving the 
HFB equation in coordinate-space can properly and self- 
consistently taken into account the deformations, pair- 
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ing, and weak-binding effects [4]. In this work, the 
Skyrme-HFB equation is solved in deformed coordinate 
spaces by the HFB-AX solver [^, where B-spline tech¬ 
niques are adopted to treat axial symmetric deformed nu¬ 
clei within a two-dimensional (2D) lattice box. The hy¬ 
brid MPI-hOpenMP parallel programming is performed 
to obtain converged results in a reasonable time [4]. 
Note that 3D coordinate-space HFB that uses multi¬ 
wavelet analysis has recently been developed [52]. For 
the particle-hole interaction channel, the UNEDFl and 
SLy4 parameterizations of the Skyrme force are adopted. 
The SLy4 force has been obtained by constraining the 
properties of neutron matter and has been widely used 
for neutron-rich nuclei [53]. The latest UNEDFl force 
is obtained by massively fitting global nuclear properties 
and is very precise for the whole nuclear landscape [53] . 
For the particle-particle channel, the density-dependent 
surface pairing interactions is employed. The pairing 
strength is taken as Uo=500 MeV for SLy4 and 460 MeV 
for UNEDFl to reproduce reasonable neutron pairing 
gaps in ^^°Sn. Note that the surface pairing has also 
been used in relativistic Hartree-Bogoliubov calculations 
for halo nuclei [25]. 

This paper is organized as follows. In Sec.2 the for¬ 
malism of the quasiparticle blocking method is intro¬ 
duced. In Sec.3 numerical calculations and discussions 
of halo nuclei ^^B, ^®B, ®^Mg are presented. Finally, the 
main conclusion of this paper and perspectives are given 
in the summary. 

2 Theoretical framework 

The HFB equation in the coordinate-space represen¬ 
tation can be written as m- 


h{r) — X 

A(r) 

' U,ir) ' 

= Ufe 

Ufe(r) 

A‘(r) 

-h{r) + X 

. . 


. . 


( 1 ) 

where h is the Hartree-Fock Hamiltonian; A is the pair¬ 
ing potential; U^. and 14 are the upper and lower compo¬ 
nents of quasi-particle wave functions, respectively; Ej^ is 
the quasi-particle energy; and A is the Fermi energy (or 

3 Results and discussion 

We have studied the ground state properties of 
and ®’^Mg with and without quasiparticle blocking. The 
density distributions, deformations, and energetic prop¬ 
erties are presented in this section. To check the de¬ 
pendence on Skyrme forces, we adopted the usually used 
SLy4 force and the newly developed UNEDFl force. 
Fig{^ displays the density profiles of ^^B calculated 


chemical potential). For bound systems, A < 0 and the 
self-consistent densities and fields are localized in space. 
For |Uj,| < —A, the eigenstates are discrete and I4(r) and 
Uk{r) decay exponentially. The quasiparticle continuum 
corresponds to \E^.\ > —A. For those states, the upper 
component of the wave function always has a scattering 
asymptotic form. By applying the box boundary condi¬ 
tion, the continuum becomes discretized and one obtains 
a finite number of continuum quasi-particles. In princi¬ 
ple, the box solution representing the continuum can be 
close to the exact solution when a sufficiently big box 
and small mesh size are adopted. In this work, the max¬ 
imum mesh spacing is 0.6 fm, the box is 27 fm and the 
order of the B-splines is 12. 

In our approach, the systems are assumed to have ax¬ 
ial and reflection symmetries. Based on the quasiparticle 
wave functions, the particle density p(r) and the pairing 
density p(r) of even-even nuclei can be written as 

P(r) = 2E.K‘(r)f4(r), 

Kr) = -‘2j:.V,{r)U:{r), 

where in the sum the quasiparticle energy cutoff is taken 
as (60 —A) MeV. With the quasiparticle blocking of the 
state fx, the density Psi^) and the pairing density Pb{'^) 
of odd-mass nuclei in the coordinate space can be written 
as [TT] : 

Pb^{v) = [/,(r)C/;(r) + U;(r)U^(r) + 2^U;(r)W(r), 

p-{v) = -2^W(r)U:(r). (3) 

Similarly, the blocking method has also been discussed 
extensively in the basis expansion method [53H55] . In 
addition to nuclear systems, the blocking method has 
also been widely used in studies of spin-polarized cold 
atomic gases [5i|3D]. The particle number equation for 
the odd-A system has to be modified as: 

N=l + 2J:k^^vl= J drps^ir). (4) 

with and without blocking method with SLy4 force and 
UNEDFl force. The densities are displayed along the 
cylindrical coordinates z-axis (the axis of symmetry) and 
r-axis (the axis perpendicular to z-axis and r=y/x'^ + y'^), 
respectively. The differences between the density profiles 
and Pr(z=o) actually reflect the surface deforma¬ 
tions. The highest state of proton lpi /2 is chosen to be 
blocked self-consistently. It shows that when calculations 
are performed without blocking, both proton and neu- 
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Fig. 1. (Color online) The density distributions of which are obtained by calculations with (a) the SLy4 force 
without blocking, (b) the SLy4 force with blocking, (c) the UNEDFl force without blocking, (d) the UNEDFl force 
with blocking. The dashed lines represent the density profiles along the cylindrical coordinates r-axis p(r, z = 0) 
and solid lines represent density profiles along the 2 :-axis p(r = 0 , 2 :), respectively. 
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Eig. 2. (Color online) Similar to Fig. 1 but for density distributions of 
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Fig. 3. (Color online) Similar to Fig.l but for density distributions of ^^Mg. 
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Fig. 4. (Color online) Similar to Fig.l but for pairing density distributions of ^^Mg obtained by calculations with 
the SLy4 force. Note that absolute values of negative pairing densities are plotted. 
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Table 1. Calculation results without or with blocking for binding energies E'“°^ and Fermi energies Xp^ and 
Xp^h, quadrupole deformations and /3“^, r.m.s. radii and pairing energies Ep°^r and Ep^f^, pairing 
gaps Ap°^ and of nuclei and ®^Mg with the SLy4 and UNEDFl forces. The subscripts wb and wob 

mean the calculations are carried out with or without blocking. Energies are given in units of MeV and radii are 
given in units of fm. 



SLy4 

UNEDFl 



19B 

37Mg 


19B 

37Mg 

j^wob 

- 98.69 

- 98.81 

- 266.14 

- 89.42 

- 90.31 

- 259.56 

^wb 

- 96.97 

- 97.20 

- 263.97 

- 87.81 

- 88.80 

- 257.62 

Qwob 

^2n 

1.80 

0.12 

^ n ( 1 . 82 ) 

2.29 

0.90 

^ n ( 2 . 37 ) 

Qwb 

^2n 

2.13 

0.23 

^„(- 0 . 58 ) 

2.44 

0.98 

^„( 0 . 43 ) 

\ wob 

Xp 

- 23.31 

- 25.32 

- 22.22 

- 22.37 

- 24.16 

- 21.46 

\ wb 
\ 

- 22.28 

- 21.79 

- 22.16 

- 21.66 

- 23.49 

- 21.17 

\ wob 

Xu 

- 1.32 

- 0.56 

- 1.98 

- 0.73 

- 0.11 

- 1.63 

\ wb 

x„ 

- 1.44 

- 0.59 

- 2.41 

- 0.79 

- 0.14 

- 1.81 

Owob 

H2 

0.001 

0.03 

0.39 

0.001 

0.001 

0.41 


0.30 

0.31 

0.45 

0.24 

0.21 

0.49 

j,wob 

P 

2.55 

2.58 

3.19 

2.52 

2.55 

3.19 

j,wb 

P 

2.55 

2.59 

3.18 

2.53 

2.56 

3.18 

j,wob 

3.27 

3.71 

3.60 

3.44 

3.80 

3.70 

j,wb 

3.25 

3.68 

3.73 

3.41 

3.78 

4.07 

j.wob 

3.09 

3.45 

3.52 

3.20 

3.51 

3.58 

j,wb 

m 

3.06 

3.43 

3.62 

3.18 

3.50 

3.84 

J^wob 

pair 

- 36.08 

- 42.50 

- 20.81 

- 31.24 

- 37.04 

- 30.55 

Tpwb 

pair 

- 28.16 

- 36.07 

- 13.47 

- 25.42 

- 31.62 

- 24.98 

A wob 
^P 

2.03 

1.89 

0.20 

1.93 

1.84 

0.01 

Awb 

^P 

10 -® 

10 -® 

10 -® 

10 -'* 

10 -** 

10 -'* 

A wob 

3.76 

3.72 

2.10 

3.42 

3.49 

2.48 

Awb 

3.61 

3.62 

1.74 

3.34 

3.42 

2.24 


tron density distributions of have a spherical shape, 
in both SLy4 and UNEDFl calculations. However, when 
calculations are performed with the blocking of proton 
lpi/ 2 , both proton and neutron density distributions of 
have obvious prolate deformation in the inner part 
and in the surface. The onset of deformation is due to 
the polarization effect caused by blocking. The total 
quadrupole deformation of ^^B with blocking is /32 = 0.30 
with SLy4. The relativistic Hartree-Bogoliubov model 
calculations of ’^^B without blocking also obtained large 
deformations [31]. In our calculations, the halo structure 
in ^^B is not so impressive. Experimentally, the den¬ 
sity distributions of ^^B can be extracted but are rather 
model dependent. The density profiles of ^®B are dis¬ 
played in Figj^ Compared to ^’^B, one can see that ^®B 
has a more pronounced halo structure with or without 
blocking calculations using SLy4 and UNEDFl param¬ 
eters. Similar to ^^B, ^®B has a prolate deformed inner 
part with blocking. Its halo is prolately deformed with 


SLy4 and is spherical with UNEDFl. 

Figj^ displays the density profiles of ^’^Mg. ®’^Mg has 
no obvious halo without blocking and shows a clear de¬ 
formed halo with blocking calculations. Recent exper¬ 
imental results demonstrated that the last bound neu¬ 
tron level of ®’^Mg has an important (about 40%) p-wave 
component [TU]. Another most likely component of the 
level is the /-wave which is close to the p-wave. A high 
centrifugal barrier of /-wave would suppress the halo for¬ 
mation. Note that in the deformed case we are blocking 
the r2’"=l/2“ state that naturally contains both compo¬ 
nents. The remarkable halo structure indicated the role 
of p-wave. We see that blocking the neutron orbit has 
little effect on proton densities of ^^Mg. Deformation of 
the inner part in ^’^Mg has little change with and without 
blocking, in contrast to 

Figj^ shows the pairing density profiles of ^’^Mg. In 
weakly-bound even-even nuclei, the pairing density dis¬ 
tributions generally have more pronounced halo struc- 
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tures than the particle density distributions. This is 
also shown in without blocking. However, the 

pairing density in ^’’Mg does not show a halo feature 
when blocking is included. Indeed, the pairing correla¬ 
tions have been significantly reduced due to the blocking 
effect. Along the z-axis, the pairing density decreases 
rapidly and goes to negative values after 17 fm, emerg¬ 
ing small amplitude oscillations. This has an analogy 
to the Larkin-Ovchinnikov (LO) phase in spin-polarized 
cold Fermi gases in elongated traps [53], although the 
spin polarization in ^^Mg is very small (~ ^). 

Table lists the calculation results of ^®B and 
^^Mg. In general, the blocking effect results in less bind¬ 
ing energies and larger deformations. It shows that 
has less obvious halo structure in calculations with SLy4 
parameters and UNEDFl parameters whether blocking 
is included or not. Experimental results [12] show the last 
two neutrons of ^^B have quite a few d-wave component 
which would hinder halo formation. The last unpaired 
proton of boron occupies the p orbit with and without 
blocking. The calculated proton r.m.s radii are almost 
unchanged in both approaches. However, the neutron 
r.m.s radii reduce slightly with blocking, as shown in Ta¬ 
ble [l] The latest experimental value of r.m.s mass radius 
of ^’^B is 3.00±0.06 fm[35]. The closest result of present 
work is 3.06 fm, from the blocking calculation with SLy4 
force. For ^®B, the experimental r.m.s radius is 3.11±0.13 
fm, which is smaller than our theoretical results. The 
two-neutron separation energy of ^®B is 0.23 MeV with 
SLy4 blocking calculations while the experimental value 
is 0.14±0.39 MeV[T3]. The total deformation increases 
and the absolute value of total binding energy reduces in 


®^Mg after the blocking. The large deformations in ^’^Mg 
are mainly due to the contributions from large spatial ex¬ 
tensions. Note that light weakly bound nuclei have very 
soft deformations. The one-neutron separation energy of 
®^Mg is 0.43 MeV with UNEDFl blocking calculations 
while the experimental value is 0 . 221 °; J 9 MeV[Tn]. One 
can see that theoretical descriptions of drip line nuclei 
are strongly dependent on the effective nuclear interac¬ 
tions. 

4 Conclusion 

In summary, the Hartree-Fock-Bogoliubov approach 
with self-consistent quasiparticle blocking in deformed 
coordinate space has been applied to study the neutron- 
rich odd-A halo nuclei ^^B, ^°B and °’^Mg. This is the first 
time such blocking calculations of odd-A nuclei in de¬ 
formed coordinate spaces has been done, and it is an ideal 
tool for describing deformed halos and continuum effects 
together. To check the dependence of the parameters, 
two Skyrme forces, SLy4 and UNEDFl, are adopted. 
Generally, the odd-A nuclei become less bound and de¬ 
formations increase due to blocking effects. In particular, 
a remarkable deformed halo structure in °^Mg has been 
obtained by blocking the U’^=l/ 2 ~ neutron state, and it 
is consistent with a large p-wave component. However, 
^^B has less obvious halos, regardless of whether or not 
blocking is used in the calculation, compared to ^°B. The 
indication of Larkin-Ovchinnikov phase in °’’Mg has also 
been discussed. 

Discussion with Xu F R, Peking University, is grate¬ 
fully acknowledged. 
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